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CN ■ ABSTRACT 
> 

. Aims. We present first results of near-infrared (NIR) J and H + K ESO-SINFONI integral field spectroscopy of the composite starburst/Seyfert 
1.8 galaxy Mrk 609. The data were taken during the science verification period of SINFONI. We aim to investigate the morphology and 
excitation conditions within the central 2 kpc. Additional Nobeyama 45 m CO(l-O) data are presented, which we used to estimate the molecular 
«j gas mass. The source was selected from a sample of SDSS/ROSAT-based, X-ray bright AGN with redshifts of 0.03 < z < 1 that are suitable 
J — for adaptive optics observations. This sample allows for a detailed study of the NIR properties of the nuclear and host environments with high 
spectral and spatial resolution. 

Methods. Integral field spectroscopy with SINFONI delivers simultaneous spatial and spectral coverage of the circumnuclear environment. 
The NIR light is influenced less by dust extinction than by optical light and is sensitive to mass- dominating stellar populations. Furthermore, 
several NIR emission lines allow us to distinguish between Seyfert and starburst activities. 
O ' Results. Our NIR data reveal a complex emission-line morphology that is possibly associated with a nuclear bar seen in the reconstructed 
^ continuum images. The detections of [Si VI] and a broad Paa component are clear indicators of the presence of an accreting super-massive 
CZ2 black hole at the center of Mrk 609. In agreement with previous observations, we find that the circumnuclear emission is not significantly 
extincted. Analysis of the high angular-resolution, molecular hydrogen emission and [Fell] emission reveals the LINER character of the 
^ , nucleus. The large H2 gas mass deduced from the CO(l-O) observation provides the fuel needed to feed the starburst and Seyfert activity in 

r> Conclusions. High angular resolution imaging spectroscopy provides an ideal tool for resolving the nuclear and starburst contributions in active 
galaxies. We show that Mrk 609 exhibits LINER features that appear to be hidden in visible/NIR spectra with larger apertures. 

Key words. Galaxies: Seyfert - Galaxies: starburst - Galaxies: fundamental parameters - Infrared: galaxies - Galaxies: individual: Mrk 609 
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1. Introduction 

Stellar kinematical studies of nearby galaxies give grounds 
speculating that every galaxy harbors a super-massive black 



hole (SIVTBH) in its center (e.g.lKormendy & Richstondll995 



Gebhardt et al. I l2000l: l Tremaine et al [2002*). Numerous studies 



of active galactic nuclei (AGN) and quiescent galaxies have 
furthermore revealed a corr elation between the SMBHs and 
the bulges hosting them (e.g. MagoiTi an et aDll998 : Page et al. 
L2OOI; Tremaine et al. 2002). These correlations indicate the ex- 
istence of feedback mechanisms that regulate the coeval growth 
of SMBHs and their associated bulges. 

In this paper we present results of an NIR study of the com- 
posite galaxy Mrk 609 (Fig.ffland IZuther et al]l2006l) with the 
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new adaptive-optics (AO) assisted integral field spectrometer 
(IPS) SINPONI at the ESO VLT. Por the first time, AO-assisted 
imaging spectroscopy on 8-10 m class telescopes allows for the 
simultaneous study of the morphology, chemical composition, 
and kinematics of the ckcumnuclear regions of AGN at an un- 
precedented spatial resolution. 

1.1. Fueling of nuclear activity 

The presence of the Seyfert phenomenon is supposed to orig- 
inate in the accretion of matter onto an SMBH in the centers 
of the galaxies. The fuel necessary for driving nuclear activity, 
which is composed of nuclear starbursts and Seyfert-like activ- 
ity, has to be transported from galactic scales (~10 kpc) down 
to nuclear scales of ~10 pc. We are still far from understand- 
ing the detailed processes that lead to the dissipation of angular 
momentum needed for the gas and stars to fall towards the nu- 
clear region. However, considerable theoretical, as well as ob- 
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servational, eff ort has been made toward understand i ng these 
proce sses (e.g. IShlosman et al. ICombeslll996t iKnapenI 
2005h . External and internal triggers of the fueling process can 
be distinguished. 

External triggers are related to the environment of galaxies 
and gravitational interaction. Non-axisymmetries, which can 
lead to loss of angular momentum, can result from galaxy in- 
teractions. For ULIRGs, which show the most extreme case of 
infrared nuclear activity, there is intriguing evidence of a con- 
nection between galaxy interaction and nuclear activity (cf. re- 



view by iKnapenI |2005| and references therein). Nevertheless, 
statistical studies show a significant fraction of nuclear active 
galaxies that are apparently free of any external trigger. Am ong 



the internal triggers, for example. IShlosman et al. (Il989h pro 



posed a two-step process that is able to sweep the interstellar 
medium (ISM) via a stellar bar from large scales into a disk of 
several hundred pc in radius. In the second step, further insta- 
bilities (bar-within-bar) drive the material close to the nucleus 
until viscous processes take over the angular momentum trans- 
port. 

While extensive star formation in quiescent galaxies ap- 
pears to be related to large-scale bars (Knapen 2005 and ref- 
erences therein), the observational evidence of bar -related nu- 
clear a ctivity is not as clear for Seyfert galaxies. i Laine et aT 



( I2OO2I) find a slight but significant increase in the galactic-bar 
fraction of active galaxies when compared to in-active galax- 
ies. This does not appear to be the case for nuclear bars. On 
the other hand, there are a considerable number of AGN that 
show no signs of any bar or inactive galaxies that do possess 
bars. The overall lack of observational evidence of direct causal 
relationships between the presence of morphological asymme- 
tries and the starburst/Seyfert activity might result from dis - 
regarding the correct spatial or time scales (i KnapenI 2005 ). 
To some extent the non-axisymmetries could occur on spatial 
scales that are currently not resolvable, or could be masked 
by dust or star formation dLaine et al ] l2002h . Considering life- 
times, the typical age for nearby spiral galaxies is 10' ° yr. Only 



10% of nearby spi r als dis play Seyfert activity (e.g. lHo et al 



1997 ). ISanders et aP ( 1988 ) propose an evolutionary scenario 
in which the Seyfert phenomenon is a transient phase. If all 
galaxies undergo this phase, the small fraction of active galax- 
ies translates into an AGN lifetime of about 10'* yr But the 
non-stellar a ctivity can also occur episodically on much shorter 
time scales ( Martini 2004). The typical bar lifetime is lO** yr 



(iCrenshaw et al.,,2003,) . an d the bar can trigger si gnificant star 
formation within 10^ yr dHunt & Malka nl ll999h . Active cir- 



cumnuclear star formation might, therefore, be related to the 
presence of a bar. But what happens with the bar when the 
Seyfert activity, under the assumption of bar- induced nuclear 
fueling, is episodic with a lifetime of about 10'* yr? Further 
work has clearly to be done, primarily aimed at the kinematics 
and dynamics of the very central region of active and inactive 
galaxies. 

The recent finding of a stellar surface mass density crite- 
rion for the separation of bulg e -domi nated and disk-dominated 



galaxies by Kauffmann et al.l (l2003h can shed some light on 
the problem. The authors find that bulges dominate in stellar 
masses ^ 3 x lO'^Mg corresponding to a mean stellar surface 



density of > 3 x 10** Mq kpd iHeckman et al.1 (|2004 find 
that this stellar surface density is also the critical surface den- 
sity for the onset of significant Seyfert activity. These results 
can be compared with the critical ste llar surface density nec- 
essary for the onset of bar instability (IWvsell2004 . It turns out 
that this surface density (y u > 1 .7 x 1 O^Mp k pc~^) is very similar 
to the one found by Hec kman et al.l (120041) for nuclear activity. 
One expects that the stellar surface density increases once a 
bar is formed, since the subsequent gas infall towards the cen- 
tral regions will be accompanied by star formation. This is seen 
in both observations and in simulations and, thus, can explain 
the factor of about 2 below the observed stellar surface density 
for bulge dominance. 

1.2. Starburst/Seyfert composite galaxies 

A s ubgroup of AGN, the starburst/Seyfert composite galax- 



ies dMoran et al. I ll996l) . appears to be best suited to studying 
the starburst-AGN connection, since the AGN and starburst 
components present themselves at the same level of activity 
(IPanessa et al .l2005h . These objects can be characterized by op- 
tical spectra that are dominated by starburst features, while the 
X-ray luminosity and its variability are typical of Seyfert galax- 
ies. The former property is based on the e missio n-line diagnos- 
tic diagrams by IVeilleux & Osterbrockl ( 1987 ). Close inspec- 
tion of the optical spectra often reveals some weak Seyfert-like 
features, e.g. [OIII] being significantly broader than all other 
narrow lines or a weak broad Ha component. Ther e is a resem- 
blanc e to narrow-line X-ray galaxies (NLXG, e.g. lBovle et al 



1995h . which also show spectra of a compo site nature. Their 
soft X-ray spectra are hard, i.e. flat (lAlmaini et al. 1996) . but it 
is still not clear how this strong and hard X-ray emission can 
be reconciled with the weak/absent optical Seyfert characteris- 
tics. The faintness of these objects in the X-ray, as well as in 
the optical domain, has not allowed us to study them in detail 
so far. 

Near-infrared (NIR) imaging spectroscopy has consider- 
able advantages over visible wavelengths: Besides the much 
smaller dust extinction, there are a number of NIR diagnos- 
tic lines (in emission and in absorption) to probe the excita- 
tion mechanisms and stellar content in galaxies. Among these 
are hydrogen recombination lines, rotational/vibrational transi- 
tions of H2, stellar features hke the CO(2-0) and CO(6-3) ab- 
sorption band heads, and forb i dden lines like JFe 11 1 and [Si VI] 
dHill et alJll999l:lMourilll994l:lMarconi etaDll994 . 



1.3. Mrk 609 

Mrk 609 shows a mixture of nuclear star formation and Seyfert 
activity and appears to be a showcase for studying the nuclear 
activity with IFSs in order to probe the coeval existence of these 
phenomena. 

There have been several observations of Mrk 609 (Figs. 
[U and [161 for the SEP). I t is a member of the original sam- 
ple of lOsterbrock with which the role of redden- 
ing in Seyfert 1.8/1.9 galaxies was studied. These observa- 
tions were followed up by X-ray, ultraviolet, optical, and in- 
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frared studies jPappa et all 120021: iRudv et al.lll988l: iGoodrich 
1990t Rudy & Rodriguez-Espinosaf 1985b . refining the clas- 



sification of Mrk 609 as a starburst/Seyfert composite. The 
UV/optical line ratios and the X-ray spectru m furthermo r e indi- 
cate a small extinction towards the nucleus. iRudv et all d 19881) 
and iGoodriciil (Il990h suggest that smaller-than-normal opti- 
cal depths and lower ionization are responsible for the large 
Balmer decrement found in the optical spectrum (Fig.fTsll. 

Optical (-0 6 yum ; iNelson etalJll996l) and NIR (-0.9 fim 



Crenshaw et al.ll2003h imaging observations of Mrk 609 with 



the Hubble Space Telescope (HST) have been used to address 
questions regarding the extension of the scattering medium 
in Seyfert 2 galaxies and the dependence of AGN fueling on 
the host morphology (Fig. [U. iNelson et al. studied a sample 
of Seyfert and non-Seyfert Markarian galaxies and found that 
Seyfert 1.8 and 1.9 show extended nuclear light profiles sim- 
ilar to Seyfert 2 galaxies. The energetics of the (circum-) nu- 
clear region of Seyfert 2s is strongly influenced by starburst 
activity. In combination with some scattering screen, this pro- 
vides the extinction responsible for the observed shape of the 
light profiles. Seyfert 1.5 galaxies, on the other hand, have 
nuclear profiles typical of Seyfert 1 galaxies, which are un- 
resolved and bright. Mrk 609 shows a strong Seyfert 1-like 
nucleus, but it appears to be slightly extended. Furthermore, 
weak broad H-recombination lines have been detected in opti- 
cal spectra , which result in an intermediate Se yfert 1.5-1.8 clas- 
sification ( Goodrich! 1990l : OsterbrocS 1981 ). This work sup- 
ports the above findings. 

Crenshaw et al. and Deo et al. I (l2006h find a larger bar frac- 



tion for narrow-line Seyfert 1 (NLSl) galaxies than for broad- 
line Seyferts (BLSl). NLSl are believed to be systems ac- 
creting close-to or above the Eddington limit. The active state 
requires an efficient fuel supply, which can be supported by 
stellar bars. As a BLSl, the large-scale optical morphology 
of Mrk 609 appears flocculent and sh ows no signs of a bar 
dCrenshaw et all 120031: lOeo et al.ll2006l) (Fig. [B- The nuclear 
morphology, however, deviate s from a circu lar symmetric ap- 
pearance. The NIR image of Nelson et alj (see inset in Fig. 
[Hi, which is less influenced by dust extinction than the opti- 
cal image, reflects this situation. Two spiral arms turn out to 
connect to an elongated nuclear structure (northeast-southwest 
direction), probably r esembling a nuclear ste l lar bar . In a re- 
cent numerical study, Englmaie r & Shlosman ( 2000 h find that 
a grand-design nuclear two-arm spiral structure can be formed 
in the presence of a weak stellar bar. Only depending on the nu- 
clear mass and the gas sound speed, such structures can allow 
for the inward mass transport beyond the inner Lind blad reso 



nance towards radii of 50-100 pc. In respect thereof iDeo et al 



(12006') find that especially the class of flocculent spirals avoids 
the presence of large-scale bars, as appears to be the case for 
Mrk 609. 

The paper is structured in the following manner. In Sect. |2] 
we describe the observations and data reduction. After a de- 
tailed analysis of the NIR and mm spectra, we discuss the 
results within this context of nuclear fueling in Sect. [3] In 
Sect. |4] we summarize the presented work and give an out- 
look. Unless otherwise noted, we use a cosmology with Hq = 




5 arcsec 



F606W 



Fig.l. HST F606W (-0.6 jum) image of Mrk 609 
The black box indicates the SINFONI 
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dMalkanet al.lll99S) 

field of view of 3 x 3 arcsec"^. The upper left inset is an 
HST F785LP (-0.9 jum) image of the nuclear 3x3 arcsec^ 
dNelson et all 1 19961) . Com pare wit h Figs. |2](d),|4](a), and the 
structure image (Fig. 2) of lPeo et al. (2006) . 



70 km s ' Mpc Q,„ = 0.3, and nA = 0.7 ( ISpergel et al 



2003h throughout this paper. 



2. Observations and data reduction 

2.1. Near-infrared data 

The observations were carried out during the science verifica- 
tion phase' (October 2004) of SINFONI, the new AO-assisted 
integral field spectrometer mounted at Yepu n, Unit Telescope 



4 of t he ESO Very Large Telescope in Chile (lEisenhauer et al 



20031) ■ The AO guiding was carried out on the nucleus of 
Mrk 609. The average seeing was around 0'.'7. The 100 mas 
pixel scale with a field-of-view (FOV) of 3"x3"was used. The 
2D image on the sky was sliced by small mirrors into 32 slitlets, 
which then were reimaged onto a pseudoslit and dispersed onto 
a 2kx2k detector. The observations covered J and H + K bands 
with integration times of 5 minutes and 30 minutes, respec- 
tively. Dispersion was achieved with the J and H + K grat- 
ings at a spectral resolution of Rj ~ 2000 and Rh+k ~ 1500. 
Successive target (T) and sky (S) observations (...TSST...) were 
acquired to produce sky-subtracted frames. 

The reduction and reconstruction of the 3D cubes were car- 
ried out using the MPE reduction software spred 3.6 ^, IDL, 
and QFitsView^. Bad pixel, cosmic rays, and flat field correc- 
tions were applied to the 2D raw frames. The 3D cubes were 
reconstructed using calibration frames for the slitlet distances 
and the light dispersion. 

Intermediate standard-star observations (near in both time 
and airmass to the target exposures) of the G2V star HIP- 



' http://www.eso.org/science/vltsv/sinfonisv/xrayagn.html 
^ Kindly provided by Matthew Horrobin (MPE, Garching). 
3 Written by Thomas Ott (MPE, Garching); 
|http://www.mpe.mpg.de/~ott/QFitsViewl 
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021070 were used to correct for strong atmospheric (telluric) 
absorptions. They can be minimized by dividing the science 
target spectrum by the standard star spectrum. This step intro- 
duces features into the resulting spectrum that are intrinsic to 
the standard star and that can be accounted for by multiplica- 
tion by the atmospheric-transmission-corrected solar spectrum. 
Since the standard star was saturated on the detector (see below 
in the context of flux calibration), we extracted the spectrum 
with telluric absorptions from the wings of the 2D image. The 
G2V spectral characteristics were removed by multiplying by 
the well known, high signal-t o-noise solar spectrum (as pro- 



vided by Maiolino et al. 199 6h. In the range of reduced atmo- 



spheric tran smission (as ind icated in Fig.|5]l, which is not cov 
ered by the iMaiolino et al.l spectrum, we interpolated the so 
lar spectrum with a blackbody, because the origi nal solar atlas 



did no t show any intrinsic features in that range dWallace et al 
1996h . 



Flux calibration is problematic, since the calibration-star 
observations saturated the detector In order to approach a 
flux calibration, igno ring the fact of variabihty (up to 3 0%) 
of Mrk 609 found bv' Rudv & Rodriguez-Espinosal(ll985 l). we 
used the 2MASS J, H, and Ks fluxes. We measured the fluxes 
within a 3"diameter aperture in the 2MASS atlas images and 
applied the fluxes to the 3" x 3" SlNFONl FOV. It has to be 
noted that this method has the drawback of ignoring the differ- 
ences in spatial resolution. We expect an accuracy of the abso- 
lute flux to be around 40%, but this is not of any significant im- 
portance, because we are primarily interested in measuring line 
ratios. For these, in principle, an absolute flux calibration is not 
necessary. We converted the magnitudes to flux densities us- 
ing the Spitzer space telescope magnitude to flux density con- 
verter'*. For 7 we measured 72mass = 13.68, which corresponds 
to//2MASS = 1.06 X 10-14 Wm-2pm-i at ^eif (^2mass) = 1-235. 
For the H + K band we used the flux density at the eff'ective 
wavelength of the Ks band. We measured Ksimk^^ - 12.50, 
which corresponds to /a:.s2mass = 4.29 x 10"'^ W m^^/im"' 
at A^s{Ks2MASs) = 2.159. The total flux within the SlNFONl 
FOV was scaled according to these values. Beca use of the low 
value of E{B - V) = 0.056 (Schlegel et al.ll 19981) . no correction 
for Galactic extinction is applied. 

From the final reconstructed 3D cube we extracted spectra 
and emission line maps. The emission line maps presented in 
Figs. |2][3] and |4] were created by summing up the flux of the 
emission line and then subtracting the average flux density of 
the continuum left and right of the line. The features seen in the 
Pao' map define the regions where we extracted spectra (Figs, 
mm 111 andO. The radius of the extraction regions is 5 pixels 
(a.'25). 

2.2. Millimeter data 

During an observing campaign of cluster galaxies with the 
Nobeyama 45m telescope in March 2005 we additionally ob- 
served Mrk 609 in '^CO(l-O). We used the autocorrelator as 
backend and an integration time of 45 minutes. The measured 
antenna temperatures were transformed to a main-beam tem- 



perature using a main-beam efficiency of 0.38. The baseline- 
subtracted CO(l-O) spectrum is presented in Fig.[T3] 

3. Results and discussion 

In this section we describe the results of the spatially resolved 
spectroscopy of the circumnuclear environment of Mrk 609 and 
the CO(l-O) observation. 

3.1. Tracing the continuum and emission line gas 

The SlNFONl imaging spectroscopy allows us to study the in- 
ner 2 kpc of Mrk 609 at a spatial resolution of about 270 pc. 
We identified several emission and absorption features in the 
individual J and H + K spectra across the FOV (see Figs. |5] 
and|7]i. Among these are hydrogen recombination lines, rota- 
tional/vibrational lines of molecular hydrogen, forbidden tran- 
sitions of [Fell] and [Si VI], as well as stellar CO(6-3) and 
CO(2-0) absorption. From selected emission lines we produced 
continuum-subtracted line maps (Figs.|2][3] and|4]l. 

In order to assess the spatial resolution we used the nuclear 
Paa emission. Fitting a 2D Gaussian to the nuclear Pao- gives 
a rough estimate for the spatial resolution, because no other 
point-source observations were available for our observations. 
We found an FWHM of about 8 pixel (~0'.'4) for the minor 
axis, while the major axis is clearly extended. We expect this 
to be a reasonable measure of the point-spread function at all 
wavelengths. This size corresponds to a linear scale of about 
270 pc (1" K 680 pc) attheredshiftofz = 0.0345. 

The morphologies of the hydrogen recombination line 
maps cleary differ from the continuum maps. The Paa map re- 
veals five emission peaks that are indicated in Fig.|2](a). The 
continuum emission presented in panel (d) of Fig. |2] shows 
an elliptical structure extending in a northeast/southwest di- 
rection and connecting to the large-scale spiral arms visible 
in the HST image (Fig. [1]). This morphology is similar to that 
of the nuclear stellar bars found in s everal other galaxies (e.g. 
Erwin & Sparic3l2002l: iMartini et al.''2003). Results of N-body 



( ,http://ssc.spitzer.caltech.edu/tools/magtojy/| l 



simulations carried out by .Patsis. (i2005i) . in which stellar par- 
ticles move on (quasi-) periodic orbits in a bar potential, also 
resemble the nuclear continuum morphology of Mrk 609. 

The Paff peaks are apparently aligned with the barred con- 
tinuum. One peak is centered on the nucleus. The four other 
peaks lie at the tip of the major axis where the bar connects to 
the spiral arms (regions 3 and 5) and on the minor axis (regions 
2 and 4). In addition to the clumpy Paa emission a fainter and 
smoother Paa component is visible, which extends between the 
distinct regions. There is a strikin g similarity between thi s con- 
figuration and the observation of lEmsellem et al.l (l2003h . who 
detected a nuclear two-arm gas spiral in the early type galaxy 
NGC 2974. This spiral structure (~ 200 pc radius) is embedded 
in a ~ 500 pc radius nuclear stellar bar. These dimensions are 
roughly the same (within a factor of 2) as for Mrk 609 (Fig. 
|2](a) and (d)). Simulations show that gaseous material can be 
transported along these spiral arms towards the central region 
(Englmaie r & Shlosm an 2000 |). Mrk 609 has als o been studied 
at 6 cm wavelengths ( Ulvestad & Wilsonlll984l) . Figure [3] (d) 
displays the Paa map overlay ed with the radio contours. The 
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radio emission appears to be slightly resolved and it also co- 
incides with the recombination line emission. Accordingly, the 
6 cm radio continuum traces the zones of ongoing star forma- 
tion. As indicated in Fig. |2](d), there is an extension towards 
regions 2 and 4, which could also be interpreted as a jet com- 
ponent. Regions 2 and 4, therefore, might be excited by the 
radio jet impinging on the inter-stellar medium. In this case the 
nuclear spiral scenario has to be altered. See Sects. 13.2. H and 
l3.5l for further discussion. 

Next, we extracted ID spectra from the H + K and J bands 
for each of the five regions by summing up the individual spec- 
tra within apertures of 5 pixel radius, centered on the Paa emis- 
sion peaks (see Figs. |5]and|2ll. The apertures sizes have been 
chosen to represent the angular resolution of the observation. 
In the resulting spectra we measured line fluxes for individual 
lines by direct measurement. Line fluxes and FWHM of the line 
profiles are given in Tab. [T] The upper part of the table corre- 
sponds to the J part of the spectrum, the lower one to the H + K 
part. Fluxes are given in units of 10"'^ W m"^ and observed 
line widths are given in km s Errors of the line fluxes were 
estimated from the noise of the neighboring continuum alone. 
One has to keep in mind that errors introduced by the flux cal- 
ibration become relevant when comparing line fluxes between 
bands (J and H + K) or when using absolute fluxes (see dis- 
cussion about calibration in Sect. 12. lb . The measurement of the 
FWHM depends on the shape of the continuum. We estimate 
its error to be on the order of 30%. In the case of Paa and Pa/3 
we fitted a broad and narrow Gaussian component with a com- 
mon line center (Figs.|6]and|8]l. The measurement of the broad 
components is quite intricate, since these lines are located in 
problematic regions of the spectrum. For example, Paff is lo- 
cated at the red end of the atmospheric transmission minimum 
( ~ 1.8;um), and the spectrum is very noisy around this fea- 
ture. As for Pa/3 and for He I, they are located at the beginning 
and at the end of the /-band spectrum and are also influenced 
by increased noise. For the nucelar l-0S(3)/[Si VI] complex, 
we fitted two Gaussians and a telluric absorption component. 
A fit without the telluric component recovers the 1-0S(3) flux 
and linewidth within the errors given in Table [T] In the case of 
the [Si VI], the flux and linewidth are considerably reduced (by 
about 60% and 30%, respectively). Thus, primarily the proper- 
ties of [Si VI] and not of 1-0S(3) are influenced. Since we only 
state the presence of the former as an indicator of nuclear ac- 
tivity and do not consider it any further, we neglect the telluric 
absorption for the 1-0S(3) measurement. 

In contrast to Ho- and HyS at visible wavelengths (FigfTsTl. 
Paa (Fig. |6]l and Pa/i (Fig. |8]l clearly show a broad compo- 
nent with an FWHM of about 3000 kms ' (after correcting for 
the instrumental spectral resolution). This supports the Seyfert 
classific ation of Mrk 609 as deduced from the X-ray prop- 
erties (iPappa et al.l l2002h . In the /-band spectrum, there is 
also evidence of a broad He I 1.083 fim component, but the 
low signal-to-noise ratio prevents us from decomposing broad 
and narrow components. Therefore, we only state the flux 
from direct measurement of the narrow component in Table [T] 
The broad component also agrees with the line width of the 
broad Lya ( ~ 360 km s ', but blended with NV) found by 



Explorer (lUE) spectral resolution of 7 A into account. The 
coronal [Si VI] emission visible in the nuclear spectrum (Figs. 
|5]and[3jc) for the line map) is another common featu re of the 



extre me energetics in AGN (Marconi et al. 1994; Pri eto et al 
I2OO5I) . 

U sing the unreddened case-B line ratio of Paa/HyS = 0.332 
dOste rbrock 1989), we find that the broad component should 
be detectable in the visible; but as mentioned in the introduc- 
tion, Mrk 609 exhibits unusual line ratios. Intrinsic extinction 
can be probed by comparing the ratios of h ydrogen recombi- 
nation lines with case-B values from lOsterbrockl (Il989h . For 
T = 10^ K, He - 10'* cm"^, and no extinction, we expect 
Paff/Bry = 12.4. Table |2] lists this ratio among others for the 
regions discussed here. Within the measurement uncertainties, 
the Paa/Bry ratios are consistent with no intrinsic extinction. 
There is in fact a trend of the ratios being a bit larger then 
the case-B value. This is especially the case for the nuclear re- 
gion, where the corresponding ratios of the broad components 
also show the same behavior. This result agrees with low H 



column densities derived from the X-ray spectra jPappa et al 



I2OO2I and t he un usual line ratios found in the UV/optical by 
Rudy et al.l (Il988h . who explained these ratios with low optical 
depths and low ionization parameters. 

The 1 .66 /im/2. 16 /jm continuum color map (Fig.fTTTi is sen- 
sitive to reddening caused by warm dust. The map shows red- 
dening along the major bar axis, with peaks at the nucleus and 
at the tip of the bar, meeting the spiral arms. Previous studies 
of the NIR continuum of Seyfert galaxies concluded that their 
J - H and H - K nuclear colors can be explained by a mixture 
of a stellar component and a warm/hot dust comp onent in emis- 
sion (e.g. H yland & Allen 1982; Glass & Mo orwoodI 1 1985 



Kotilainen et al.lll992l; lAlonso-Herrero et al.lll996 ). Using flie 
H and K zero points from Table 3.4 in lGlasslll999l) . \hQH-K 



colors range from 0.6 in the outer regions up to 1.0 at the 
nucleus. These are compatible w ith warm (~ 500 K) dus t 



em ission and small ext incti on (cf. [ Glass & Moorwoodlll985h . 
As lCalzetti et al.1 (Il996l) and lFischera et al.l (l2003b furthermore 
show, the effects of dust in starbursts can often be fairly well 
modelled by a homogeneous or clumpy/turbulent for eground 
screen. Since the extinction curve is very flat in the NIR (IDraine 



19891) and since the wavelengths of the lines used in the present 



Rudy et al.l (119881) . when taking the International Ultraviolet 



analysis have small separations, the differential reddening be- 
tween these lines will be small. Therefore, the H-K reddening 
will only have a neglegible influence on the measured line ra- 
tios. 

Due to the uncertainties in the absolute calibration, where 
we do not have reliable knowledge of the host galaxy and nu- 
clear contribution in the measuring aperture, we will not use 
the Pa/)/Pa/3 ratio as an indicator for extinction. Also the J 
band spectrum shows no significant extinction. If we assume 
the same case-B conditions as above, we would expect not to 
detect the Pay line (close to He I 1.0883 fim). Figure |4] how- 
ever, shows faint signs of Pay, which points towards the same 
conclusion of small extinction around the nucleus. 

Pap, Bry, and He I 2.058 /im exhibit the same spatial dis- 
tribution (Fig|2](b), (c), and Fig.|4](d)), and their intrinsic line 
widths are comparable and are spectroscopically unresolved. 
This indicates that the lines arise in the same parcels of gas. 
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E 




Fig. 5. H + K rest-frame spectrum extracted from a 5-pixel ra- 
dius aperture centered on the nucleus (region 1 in Fig.|2](a)). 
Prominent lines are indicated (cf. Table [T]). The high noise 
around 1.8 fim, indicated by the thick bar, is due to reduced 
atmospheric transmission between the H and the K bands. 



E 




1.10 



1.15 



1.20 



1.25 



1.30 



K., [Mm] 



Fig. 7. J rest-frame spectrum extracted from a 5 pixel ra- 
dius aperture centered on the nucleus (region 1 in Fig.|2](a)). 
Prominent lines are indicated (cf. Table[T]). 



1 6.0 




Fig. 6. Rest-frame spectrum around Paa, extracted from the nu- 
clear 5 pixel radius aperture. The spectrum has been smoothed 
with a 3 -pixel boxcar The blue dashed and red solid lines rep- 
resent Gaussian fits to the narrow and broad components of 
Paa. 



Note that on the nucleus Bry is only measured as an upper limit 
and He I 2.058 jim only has a narrow component. On the other 
hand He I 1.0883 yum is just concentrated on the nucleus (Fig. 
|4](b)) and has both a broad component and a narrow one. 



3.2. Molecular hydrogen and iron 

The H2 emission lines originate in surfaces of molecular clouds 
exposed to stellar or nuclear radiation. The observed line 
strengths and ratios strongly depend on the excitation mech- 
anisms of H2 discussed below. 

[Fe II] emission is believed to originate in partially ionized 
zones, which occur around supernovae or the active nucleus. 
Such regions can be produced by the hard ionizing X-ray/UV 



E 
a. 




Fig. 8. Rest frame spectrum around Pa^, extracted from 
the central 5 pixel radius aperture. The spectrum has been 
smoothed with a 3 -pixel boxcar and the result of the fit of a nar- 
row (blue dashed curve) and a broad (red solid curve) Gaussian 
component with a common line center is shown. 



Table 2. Ratios of prominent emission lines using the narrow 
components. 



Region 


Pair 
Bry 


H2i>.li!l,im 
Bry 


|t'ellJl.i!i'Vm 
Pa/3 


1 


> 14.3 (> 13.7)" 


> 17.3 


3± 1 


2 


15 ±3 


3.2 ±0.7 


0.8 ±0.2 


3 


12 ±2 


1.0 ±0.2 


0.32 ± 0.04 


4 


14 ± 3 


2.3 ±0.6 


0.5 ±0.1 


5 


12.2 ±0.1 


0.6 ± 0.2 





Broad component in brackets. 



continuum or by sh ocks in the circumnuclear g as or in super- 
nova remnants (e.g. lAlonso-Herrero et al ] |l997h . 
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Table 1. Emission line fluxes for the 5 regions identified in the Paa line map (Fig.|2](a)). 



Line 


Region 


1 


Reg 


ion 2 


Region 3 




Region 4 


Region 5 




A [p.m] 


Flux 


FWHM 


Flux 


FWHM 


Flux FWHM 


Flux 


FWHM 


Flux FWHM 


He I /II. 0833 narrow 


31.8 ± 10.6 


540 


— 


— 


— 


— 


— 


— 


— 


— 


[Fell] i 1.257 


5.1 ± 1.0 


439 


1.3 ±0.3 


327 


0.7 ±0.1 


160 


1.3 ±0.2 


303 






Pa/S narrow'' 


1.9 ±0.6 


214 


1.7 ±0.1 


150 


2.2 ±0.1 


163 


2.2 ± 0.2 


183 






PayS broad" 


15.0 ±7.0 


3400 


















[Fell] AIM4'= 


2.7 ±0.2 


410 


0.6 ±0.1 


240 


0.3 ±0.1 


200 


0.5 ±0.1 


260 


0.07 ± 0.01 


200 


Paa narrow 


5.3 ± 0.6 


300 


7.9 ± 0.2 


270 


7.0± 0.1 


200 


5.4 ± 0.2 


200 


0.22 ± 0.01 


210 


Paff broad 


16.4 ±5.6 


3000 


















l-0S(3)''-= 


5.0 ±0.8 


317 


1.1 ±0.2 


371 


0.5 ±0.1 


345 


0.6 ± 0.2 


471 


0.05 ± 0.02 


230 


[SiVI]'''^ 


16.4 ±4.2 


2130 


















1-0S(2) 


2.1 ±0.4 


308 


0.6 ±0.1 


366 


0.20 ± 0.06 


321 


0.4 ±0.1 


333 


0.06 ± 0.02 


450 


He 1^2.058 


0.3'' 


300 


0.3 ±0.1 


338 


0.20 ± 0.05 


200 


0.20 ± 0.05 


200 


0.07 ± 0.01 


200 


2-18(3)= 


0.2 


317 


0.1 


322 


0.1 


312 


0.1 


316 


0.02 


300 


l-OS(l) 


6.4 ± 0.4 


330 


1.6 ±0.1 


375 


0.6 ± 0.04 


270 


0.9 ±0.1 


300 


0.10 ±0.01 


300 


Bry 


0.37" 


300 


0.5 ±0.1 


267 


0.6 ±0.1 


200 


0.4 ±0.1 


200 


0.18 ±0.01 


200 



" Decomposition using two Gaussian components. 

1-0S(3) and [Si VI] are decomposed using two Gaussian components. 

Result strongly influenced by telluric absorption. 
" Upper limit, assuming the FWHM of the narrow Pao- component. 
° Upper limit, assuming the average FWHM of the other molecular hydrogen lines. 



3.2.1 . Kinematics of tine H2 and [Fe II] lines 

The H + K spectra exhibit a couple of rotational/vibrational 
lines of molecular hydrogen (Fig. |5]l. The most prominent is 
the l-OS(l) transition. Figure [Sfb) shows the spatial distribu- 
tion of H2, which appears rather confined to the nucleus. Most 
of the l-OS(l) emission is concentrated within 260 pc (30% 
peak contour in Fig. lUa)) around the nucleus. The emission 
further shows a smooth decline towards the outer regions as 
it follows the continuum contours. This is in agreement with 
previous investigations of the molecular content of samples of 
Seyfert galaxies, which showed that the gas is concentrated in 



a disk around the centers of th e galaxies dReunanen et al .I2OO3I 



Rodriguez- Ardila et al.ll2005l) . The [Fell] emission is oriented 
almost perpendicular to the H2 emission and follows the recom- 
bination line gas (Figs.|3](b) and|2](a)-(c)). Comparison of the 
line width of H2 with the one of the naiTow lines can thus yield 
constraints on the location of the molecular and narrow-line 
gas. We consider line widths to be spectroscopically-resolved, 
if the FWHM is larger than 280 km s"', i.e. V2 times the spec- 
tral resolution of 200 km s '. From Table [T] we find that the 
molecular hydrogen emission is slightly resolved spectroscop- 
ically and that the line widths are similar in all regions. The H2 
lines are also always broader than the narrow recombination 
lines. In contrast to this, [Fe II] 1 .257 fim appears to be resolved 
spectroscopically only on the nucleus and in regions 2 and 4. 
In these three regions [Fe II] is also broader than the narrow re- 
combination lines and even broader than molecular hydrogen. 
Apart from that its FHWM is similar to that of the narrow re- 
combination lines. Off-nuclear [Fell] follows the patchy Pao- 
emission/Bry/He I distribution, thereby tracing star formation. 
On the other hand, the unusual broadening of the iron emission 
in regions 1, 2, and 4 might also be related to a jet-like exten- 
sion (southeast-northwest direction) in the 6 cm VLA contin- 



uum map of lUlvestad & WilsonI (Il984 . Figure [3] (d) displays 
our Pao' map overlayed with the radio map in contours. It is 
rather speculative but conceivable that a radio jet might alter the 
iron gas kinematics (cf. discussion in Reunanen et al. 2003)). 
These observation s se em to confirm previous result s by 



Reunanen et al 



2003^) andll odrfguez-Ardila et all (12005!) that 



the H2 kinematics is decoupled from the [Fe II] and narrow re- 
combination line kinematics. 

However, because of the combination of low inclination 
and limited spectral resolution, we cannot detect any consid- 
erable dynamics from line shifts. 

3.2.2. H2 excitation mechanism by line ratios. 

Three distinct processes are commonly discussed: (1) UV flu- 
orescence (non-thermal), (2) X-ray heating (thermal), and (3) 
shocks (thermal). The processes result in different H2 responses 
and can thu s helping dis tinguishing the dominant excitation 
process (e.g. lMourilll994 . 

(1) Excitation via UV-fluorescence can occur via UV pump- 
ing and/or collision with H-atoms, which are controlled by tem- 
perature, density, and strength of the UV radiation field. Within 
a warm, high-density gas in a strong UV radiation field, ther- 
mal line ratios are found for the lower vibrational levels. In this 
case the lower levels are dominated by collisional excitation, 
whereas the higher levels are populated via UV pumping. To 
distinguish this scenario from pure shock excitation, observa- 
tions of the high transitions are essential. In the UV excited gas, 
strong 2-IH2 lines are expected, as well as strong lines in the H 
band, e.g. 6-4Q(l) 1.6 um o r 1-0S(7) 1.75 urn (cf. Table 2 of 



Black & van Dishoeck 



19871) . The //-band lines only become 



weak relative to the /T-band lines in the case of very high den- 
sities and high UV intensities. Since we detected no H2 lines 
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in the H band and the 2-IH2 lines are only measured as upper 
limits, UV radiation as the main excitation mechanism seems 
unlikely. 



(2) Excitation by X-r ays or cosmic rays. iTine et al 



Draine & Woods (Il990l) . as well as iMalonev et al 



1997) 



1996) 



analyzed the excitation of H2 by X-rays and cosmic rays. 
Obvious non-thermal line spectra from these X-ray-dominated 
regions (XDRs) are only observable, if the temperature is 
well below 1000 K and the ionization-rate per H-atom 
is below 10^'^cm^s"'. At higher temperatures, the lower 
vibrational levels will be populated mainly by collisions, 
whereas higher ionization-rates will destroy the H2 molecules. 
iLepp & McCray ( 1983) derived the relation L/r x lO^^Lx be- 
tween the 1 - 10 keV X-ray luminosity and the l-OS(l) line lumi- 
nosity for isobaric models containing a compact X-ray source. 
The 2-10 keV X-ray flux of Mrk 609 (Pappaet al. 2002) is 
about 2.7 X 10"'^ W roT^ and translates to a l-OS(l) line flux 
of about ~ 10"'^ W m This is comparable in magnitude to 
the measured l-OS(l) flux, if summing over all regions. The 
analysis of the line ratios, however, does not strongly support 
excitation via X-rays. 

(3) Thermal excitation via shock fronts. Thermal excitation 
occurs through collisions with H or H2. The population of the 
electronic ground levels in the ro-vibrational transitions repre- 
sents a Boltzmann distribution. Then the temperature of the gas 
represents a kinetic temp erature that could be as hi gh as 2000 K 
or more. See for example Draine & McKee ( 1993^) for a review. 

In order to test whether the [Fe II] and H2 emission can be 



caused by X-ray heating, we used the models of Malonev et al 
(11996') to estimate the emergent l-OS(l) and [Fell] fluxes of 
a gas cloud illuminated by a source of hard X-rays with an 
intrinsic luminosity Lx. The cloud is at a distance d from the 
X-ray source and has an electron density n^. The cooling is then 
given by the effective ionization parameter ^eff : 



^eif = 1.26X 10" 



fx 



n,N°:' 



(1) 



'22 



where fx is the incident hard X-ray flux at the distance d [pc] 
from the X-ray source, ns [lO^^cm^^] is the total hydrogen gas 
density, and A^22 [lO^^cm^] the AGN intrinsic attenuating col- 
umn density. The va lue of fx can be calculated via Lx/4-7rd^. 
Maloney et al. ( 19961) calculate the emergent intensities for two 



gas densities, 10 cm" and 10 cm" , which can be directly 
read off their Fig. 6(a, b). Then, the flux can be calculated us- 
ing the solid angle of the cloud region, i.e. the aperture of a 
10-pixel diameter (O'.'S corresponding to 4.6 x 10"'^ sr). The 
BeppoSAX 0.1-10 keV X-ray luminosity is 6.3 x 10"*^ erg s"' 
(in the high state) and the absorbing column den sity is Nh < 
1.32 x 10^' cm"^ for a single power-law model (Papp a et al 



20021). For the indicidual regions, Table[3]lists the H2 and [Fell] 
fluxes in units of 10"^" W m"^ for the distances li in pc. 

Taking the uncertainties in our flux calibration into account, 
the H2 estimates show that X-ray heating can account for part 
of the observed l-OS(l) flux (cf. Table [U in the nuclear re- 
gion. The model also predicts the 2-lS(l) and the Bry fluxes. 
The expected line ratios approximately resemble those of the 
observed 2-lS(l)/l-0S(l) and 1 -OS (IVBry ratios. On the other 



Table 3. Emergent H 7 and [Fe II] fluxes using the models of 
IVIalonev et al. I (ll996l) . 



d 


n = 10^ 
l-OS(l) 


cm ' 
[Fell] 


n = 10' 
l-OS(l) 


cm ' 
[Fell] 


50 


1450.0 


0.5 






80 


15.0 


0.0 


460.0 


145.0 


420 


15.0 


0.0 


0.05 


0.0 


630 


9.0 


0.0 


0.08 


0.0 


1300 


3.0 


0.0 


0.15 


0.0 



hand, [Fe II] 1 .644 fim appears not to be excited by X-ray heat- 
ing, except for the nuclear region in case of the lower electron 
density environment. This indicates that other processes like 
shock heating or fluorescence have to be accounted for. In fact, 
the NIR diagnostic diagram (Fig.[T4b points towards the impor- 
tance of shock heating. 

In summary, the observed line ratios displayed in Fig. |9] 
agree quite well wi th those of a shock model with 2000 K from 
Brand et al. (Il989l) for all regions. A scenario with mixed ex- 
citation is possible and more realistic as well. Fitting the 2-1 
S(l)/l-0 S(l) ratio with the models 14 (UV excitat ion) and S2 
(shock excitation) of lBlack & van Dishoeckl (Il987h for regions 
1-4 about 96% of the line emission should be due to thermal 
excitation and about 4% due to UV pumping. For region 5 ther- 
mal excitation could account for up to 60% and UV pumping 
for up to 40%. 



1.5 



0.0 




3000 K 



_i_ 



non-thermal 



_i_ 



0.0 



0.2 0.4 
2-1S(1)/1-0S(1) 



0.6 



Fig. 9. H2 diagnostics diagram using the two line ratios 
2-lS(l)/l-0S(l) vs. 1-0S(3)/1-0S(1). The region with the 
horizontal fill pattern covers the non-thermal models of 



Black & van Dishoeckl (119871) . The vert ically patterned region 
are the thermally uv-excited models of [Sternberg & Dalgarnd 
19891) . T he open triangle corresponds to the X-ray heating 



shock model of Brand et al. 



models of iDraine & Woodsl (11990) and the open circle to the 



1989). Filled squares are the mea- 
sured line ratios with numbers indicating the extraction region 
(see Fig.|2](a)). 
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3.2.3. The H2 population diagram. 

Another approach to finding the relevant excitation mechanism 
is to assume thermal excitation. The ro-vibrational levels are 
then populated according to the Boltzmann equation. This can 
be illustrated by a population diagram showing the observed 
population density versus the energy of the upper level. If all 
values lie on a straight line, the gradient is proportional to 
the temperature of the molecular gas. The population density 
can be derived from the observed column density according to 
Lester etalJ(ll988h : 



f 

N,o\ = — — X 



An 



Aui he Q; 



(2) 



aperture 



where [/] = W m ^ is the observed line flux, Aui the transi- 
tion probability (taken from lWolniewicz et al ]ll998h . A the rest 
frame line wavelength, h the Planck constant, and c the velocity 
of light. The additional factor 47r/Qapeiture takes the spatial dis- 
tribution into account. The Boltzmann distribution for thermal 
equilibrium can be written as 



N" g'J 



(3) 



where A^' and A^" denote the column densities at the corre- 
sponding levels, g'j and g'J are the statistical weights, and T 
the temperature of the thermal equilibrium. If the H2 is non- 
thermally excited, only lines within one rotational level (e.g. 
V = 1 - 0) fall on a straight line. 




region 1 

region 2 

region 3 

region 4 

region 5 
-linear fit to region 1 data 
_] I I I I I I 



v = 2-1 



_i_ 



6000 7000 8000 9000 10000 11000 12000 13000 14000 

excitation temperature [K] 

Fig. 10. H2 level population diagram relative to l-OS(l) for the 
the five regions of Fig. |2] (a). All v = 2 - 1 fluxes are upper 
limits. A linear fit to the region 1 data is presented and yields 
an excitation temperature of about Tex - (1631 + 107) K (cf. 
Fig.©. 



The population diagram Fig.[TO]for the H2 lines of Mrk 609 
shows that, within their errors, the different transitions fall on 
a straight line. From this we estimate a kinetic temperature for 
the nuclear region of about 163 1 + 107 K. It has to be noted that 



the V - 2 - I fluxes are upper limits. For region 1 these values 
are consistent with the v - 1-0 kinetic temperature. For re- 
gions 2, 3, and 4, the scatter with respect to a linear relationship 
is larger and the v - 2 — \ values tend to show lower excita- 
tion temperatures. Nevertheless, the temperatures found in the 
population diagram agree with the temperatures indicated by 
the solid curve in Fig. |9] A stronger non-thermal contribution 
might be possible for region 5. 




Fig. 11. Color map derived from continuum fluxes measured 
at the 2MASS eff'ective H and Ks band wavelengths and zero 
points. Darker greyscales indicate redder colors with H - K 
ranging from 0.5 to 1. Contours correspond to Pao- emission 
displayed in Fig.|2](a). 



3.3. Stellar absorption features 

The H + K spectra of all five regions show stellar absorption 
features like the CO(6-3) 1.62 jum and CO(2-0) 2.29 jum, as 
well as Si I 1.59 /im, Mgl 1.5,1.71 /im, and NaI2.206//m (see 
Fig.fTSI). The details of the stellar content will be discussed in 
a subsequent paper, but one can get initial insight into the stel- 
lar populations by using the ratio of equivalent widths (EW) 
of CO(6-3) an d Si I 1 .59 urn to i ndicate the t emperature of 
late-type stars dOridia et al]|l993l: llvanov et al]r2004) . CO(6- 
3) grows rapidly from early K to late M stars, while Sil varies 
only very slowly with temperature. Table|4]lists the ratio for all 
regions. The measured ratios co rrespond to M-type giants (see 
Fig. 5(b) in lOriglia et al.lll993h . One has to take into account 
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that in active galaxies these stellar features can be substantially 
diluted by non-stellar nuclear emission. This effect can be rec- 
ognized by comparing spectra of regions 1 and 2 of both panels 
in Fig. [12] While being very similar in the H band, the nuclear 
/T-band spectrum is significantly reddened. However, we do not 
go on to consider the effect at this point. 



Table 4. Measurements of stellar absorption Unes. 




1.60 



1 7 \ 



INqI 


lC0('2-0) 








.-~---vvv^^.Ai\ region 1 






region 4 


'■'■'''■'■''"''^^^ region 3 


- 1 region 5 





2.22 2.24 2.26 2.28 2.30 



2.32 2.34 



Fig. 12. Stellar absorption features in the //-band {upper panel) 
and /T-band {lower panel) parts of the spectra extracted from 
the five regions of Fig.|2](a). The fluxes of the H and K spectra 
are scaled to the flux at 2.1 yum of region 1. Then the spectra 
were shifted by an arbitrary amount for better visibility. The 
/T-band part of region 1 between 2.225 and 2.275 yum is not 
useable because of a detector defect and has been replaced by a 
flat line. The same situation applies for region 5 between 2.241 
and 2.27 //m. 



The CO(2-0) inde x can be used to esti mate the age of 
the stellar population ( Origlia & Olivall200d) . As outlined in 



Ramos Almeida et al 



( 200d) . 



we use the Ivanov CO index 
(llvanov et al., ,20041) and convert it to a photo metric CO index 



to rea d the approximate age from Fig. 1 of lOriglia & Oliva 
(l2000l) . The photometric CO indices range from 0.1 to 0.15. 
For solar metallicity this corresponds to ages between 80 and 
150 Myr for the Geneva tracks with the younger ages mea- 



Region 


EW(CO(6-3))/EW(SiI) 


CO(2-0) index 


EW(Nal) 
[A] 


1 


0.31 


0.09 


2.5 


2 


0.28 


0.19 


2.8 


3 


0.26 


0.14 


3.0 


4 


0.31 


0.13 


3.0 


5 


0.24 


0.17 


3.8 



sured in regions 2 and 5. Such d eep CO absorption i s typically 



sor ption 1 

found in starburst galaxies (e.g. iFischer et alJ [ l2006l and refer- 
ences therein). Note that the blue edge of the CO(2-0) bandhead 
is influenced by atmospheric OH emission, and the change in 
index is on the order of 20%. Thus, the stellar ages might be- 
come shorter by a few 10 Myr 

The Na 1 absorpti on is commonly fou nd in late-type stars. 
Model calculations bv lDavies et al.l ( l2005b show that the EW of 
Na 1 typically ranges between 2 and 3 A. Their models include 
AGB phases that have a si gnificant influence on the depth of the 
Na I feature (see Fig. 7 of iDavies et al ]l2005h . Indeed, the EWs 
range between 2.5 (for region 1) and 3.8 (for region 5), which 
is consistent with the ages found from the C O(2-0) index. Both 
results agree with the young stars found by iRudv et al.l (119881) 
from their UV/optical spectroscopy. 

3.4. 00(1-0} observations 

Another ingredient in the starburst/Seyfert picture is the pres- 
ence of fuel in the form of molecular gas. Giant molecular 
clouds are the birthplaces of stars. The dense gas is dissipa- 
tive and galaxy interaction is believed to channel huge amounts 
of molecular gas towards the nucleus of the int eracting galax 



ies, p roviding fuel for the nuclear activity (e.g. lSpringel et al 



20051) . Infrared luminous AGN { Ljr > 10 " erg s ') are fo und 



2001 



20051 and 

10i°M„. 



to be rich in molecular gas (e.g. lEvans et al 
references therein) w ith gas masses of up to M^^ 
Maiolino et alj (1 1997b explain the more intense star formation 
found in Seyfert 2 galaxies with a significantly higher frac- 
tion of asymmetric morphologies than in Seyfert 1 and field 
galaxies. As discussed in Sect. 11.11 there is compelling ob- 
servational and theoretical evidence that bars efficiently redis- 
tribute angular momentum in galaxies and cause gas to flow in- 
wards into the circumnuclear reg ion (inner 1-2 kpc), giving rise 
to starburst and n uclear activity dKormendv & Kennicuttll2004l: 
Joge e et al.ll2005b . A major result from the Nuclei of Galaxies 
(NUGA) project carried out at the Plateau de Bure interferom- 
eter is the wide variety of circumnuclear d isk morphologies on 



the 10 pc scale ass ociated with the A GN dGarcfa-Burillo et al 
20051 and also ,Schinnerer et al.ll2000l for an integrated CO/NIR 
study of NGC 1068). It is believed that on these scales, sec- 
ondary perturbations modes appear to take over and are respon- 
sible for channeling gas towards the center of the galaxy. Still 
the scales are large compared with nuclear scales. The actual 
fueling of the central engine is believed to be supported by vis- 
cous flow of the material delivere d by the secondary perturba- 
tions dOuschl & Strittmattej|2004l) . 
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Fig. 13. Observed-frame CO(l-O) spectrum of Mrk 609. 



During an observing campaign of cluster galaxies with the 
Nobeyama 45m telescope in March 2005, we also observed 
Mrk 609 in '^CO(l-O). The baseline-subtracted, main-beam- 
efficiency corrected CO(l-O) spectrum is presented in Fig. [T3] 
The beam with a size of about 15"essentially covers the entire 
visual part of the galaxy (cf. Fig. [B. The FWHM of the CO 
line is 51 +4 km s ' and the line-integrated flux, corrected for 
main-beam efficiency, is Iqo = (14 ± 1) K km s"'. The narrow 
line width is reminiscent of the small inclination of Mrk 609. 



According to Eq. 2 in lSolomon et al.l (11992) . the line luminos 



ity results in L^q = (1.5 + 0.1) x lO' Kkm s"> pc^. We can now 
estimate the H2 mass, assuming optically thick and thermal- 
ized emiss ion, originating fro m gravitationally bound molecu- 
lar clouds (IE vans et al ] l2005h : 



M(H2) 
^co 



(Kkms^'pcV 



(4) 



The c onversion factor a lies between 2 and 5 (iRadford et al 
199 lb . Here we adopt a - 4, resulting in M(H2) = (6.0 + 0.4) x 
lO'^Mo, which is on the high end of t he range of masses found 
in AGN (e.g. lRigopoulou et al.|[l997l) . 

The gas mass can be compared with the dust mass derived 
from the IRAS 6 fim and 100 um fluxes (fm ^ 2 550 mJy, 
/loo = 4760 mJv; lRudv & Rodriguez-Espinosalll985i) . The in- 
frared emission originates from warm dust, heat ed either by 



star-fo rmation or by the active nucleus. Following lEvans et al 
(l2005b . the dust temperature results in Tdust ~ 47 K, and 
the corresponding dust mass in Mdust ~ 1.1 x lO^M©. The 
gas-to-dust ratio, ~ 545, is compar able to /jM5-detected spi - 



ral and luminous infrared galaxies (lYoung & Scovilldll991l) . 
Next, w e calculate the infrare d luminosity using the formula 
given in lSolomon et al. ( 1997 ). Tutui et al. ( 2000l) furthermore 
show that an IRAS /T-correction for this type of galaxy is not 
necessary. The color correction, r, lies between 1.5 and 2.1 
dSolomon et al .Il997h . Assuming r - 1.8, the infrared luminos- 
ity calculates as Lfm lO^L©, which is typical for an infrared 
luminous galaxy. 



The FIR luminosity of (inactive) galaxies is interpreted as a 
measure of the number of visible and UV photons, thus measur- 
ing the number of high-mass stars. This allows us to estimate 
the global star-formation rate (SFR) from the IRAS fluxes us- 
ing Eqs. 7 and 8 in lHopkins et all (l2003l) . The global SFR of 
Mrk 609 amounts to ~ 30Mq yr"', on the other hand, CO(l- 
0) measures the cold molecular gas supplying the star forma- 
tion. The ratio between FIR and CO luminosity, t herefore, is a 



measure of the star formation efficiency (SFE; cf. lEckart et al 



1994). The global SFE is on the order Lfir/Mh, * 17 Lq Mg\ 
which is again on the high end of the SFEs of non-interacting 
galaxies. Suc h SFRs and SFEs are als o often found in interact- 
ing systems (iSolomon & Sagelll988l) . Stars thus appear to be 
formed very efficiently in Mrk 609 and the gas depletion time 
is short (~ 3 X 10** yr). 

We can also estimate an SFR using the 21 cm radio flux. 
For Mrk 609, the NRAO VLA Sky Survey (NVSS) 21 cm 
flux /21cm - 30 mJy yields a luminosity of Lgicm - 8.3 x 
10^^ W Hz According to Eqs. 1 and 2 of iHopkins et al 



(l2003h . the radio derived SFR amounts to ~ 46 M[7) yr . This 



is som ewhat higher than the FIR-derived SFR. iHopkins et al 



(12003b find a tight correlation between radio and FIR flux. The 
differences in the derived SFRs could be caused by a stronger 
AGN contribution (e.g. jet) at the 21 cm wavelengths. 



3.5. The starburst/AGN connection in Mrk 609 

In contrast to previous observations of Mrk 609, we have now 
spatially resolved the inner circumnuclear environment on the 
270 pc scale. This allows us to study the relative contributions 
and importance of nuclear and off-nuclear (HII regions, late- 
type stars) emission. 



Analogously to the optical (e.g. 'Baldwin et alj 1981 : 
^illeux & Osterbrock 1987), Rodriguez-Ardila et all J2004 . 
20051) emphasize the use of the line ratios [FeII]/Paj8 and 
1-0S(1)/Bry to distinguish between starbursts, AGN, and 
LINERs (low ionization nuclear emission-line region galaxy). 
They find that objects with either one of the ratios being lower 
than 2 were predominantly identified as Seyferts. Starburst 
galaxies, on the other hand, are located in a region where both 
ratios are lower then 0.4, while LINERs have both ratios higher 
than 2 (Fig. [Hi. 

In this context, the ratios describe the transition from pure 
shock excitation driven by supernova remnants (upper right 
corner of the diagram) to purely ionizing radiation powered 
by star formation (low er left corne r). For illustration purposes 
we plot the data from iRodriguez-A rdila et al. (2005) as open 
symbols in Fig. [14] Starbursts from iPale et al.l (120041) populate 
the l ower left region of the diagram . Seyfert 1 and 2 galaxies 
from iRodrfguez-Ardila et al. j 2005 [populate the middle part. 



whereas LINERs dLarkin et al.lll998h can be found in the up 



per right corner. The emission of the supernova IC 443, on the 



other hand, is believed to be driven by shocks (IGraham et al 



19871) . The location of AGN (powered by X-ray heating) be- 
tween these two extremes has not been fully understood yet. 
But most probably it is the result of a power-law continuum il- 
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Fig. 14. Line ratios of [Fell] 1.257;um/Pa^ and 1- 
0S(1) 2.121//m/Br7. Activity types (starburst, AGN, LINER) 
are indicated. Filled symbols represent our measurements 
of the five regions of Mrk 609, as well as for the total 
FOV (red star). Because of missing /-band data, region 
5 is only accurately located along H2/BTy. Open symbols 
correspond to hterature v a lues: rhombs represent starburst 
galaxies from Dale et al. (200 41), left tria ngles Seyfert 1 
galaxies from Rodri g uez- Ardila et alj (|2004. headlong trian- 



gles Seyfert 2 s from iRodriguez-Ardila et al.l (I2005h . circles 



LINERs from "Larkin et alJ 
also from Larkin et al.. (199 



19981) . and asterixs supernovae 



. The diagram is interpreted as 
displaying the transition from pure photoionization (lower 
left corner) to pu re shock driven (upper right) emission 
dRodriguez-Ardila et al. 2005). Assuming a linear relationship 
between the two line ratios, the dotted magenta line represents 
a linear fit to the literature data points. The blue dashed lines 
represent the Icr prediction band of the fit. 



luminating a slab of gas (see below and iRodrfguez-Ardila et al 
2005h. 



It can be seen from Figs. |2] (b) and |4] (c) that the [Fell] 
emission is extended towards regions 2, 3, and 4. This extent 
is about perpendicular to the H2 one. Moreover, the radio con- 
tours are approximately congruent with the [Fe II] distribution. 
The NIR line-ratios measured in the five regions fall into dif- 
ferent parts of the NIR diagnostics diagram (Fig.[T4li. The cali- 
bration uncertainties discussed above play no role at this point, 
since the line ratios are calculated within a spectral band. Our 
data points clearly follow a trend from a Liner-like value at 
the nucleus to a starburst-like value at the most distant region 
5. The circumnuclear regions 2, 3, and 4 fall in the domain 
of potentially mixed excitation. We discuss the situation for 
the individual regions below. Under the assumption that the 
literature values follow a linear relationship between the line 
ratios (from photo ionization- to shock-driven excitation, see 
Larkin et"al]|l998h . we carried out a linear regression shown 
in Fig. [T4l (correlation coefficient R ~ 0.6 and a probability of 
= of about 10 '^) with its Icr prediction band. This fit, how- 



ever, is only intended to guide the eye. A firm correlation has 
not been established yet. Notice that our data fall somewhat 
below the linear relation, but still within the prediction band. 
Incorrect estimations of the continuum levels during the mea- 
surement of emission line fluxes might in part be responsible 
for this trend. On the other hand, M-type giants are the dom- 
inating population in the NIR (Sect. 13.3b . and stellar absorp- 
tion can t herefore influence the flux measurement of emission 
lines (e.g. i Harrison et al. 1998 1 who found that stellar absorp- 
tion play s an importa nt role in their observation). The study 
of iLarkin et al.l (Il998h used larger apertures compared to our 
work from which they extracted spectra. In this case a consider- 
able part of the host galaxy emission c an modify the emission 



hne measure ments. On the other hand. [Rodriguez- Ardila et al 



( I2005L 120041) use nuclear apertures of sizes comparable to our 
work. The latter study, however, does not consider circumnu- 
clear stellar absorption features at the position of the emission 
hnes, which can influence the results. 

Visual inspection of the spectr a of late-type stars 
dWaUace et al.ll2000t llvanov et al1l2004l in J and H, K, respec- 
tively) shows that there is no significant stellar absorption at 
the position of the emission lines. Atmospheric residuals might 
account for an uncertainty in flux on the order of 20%, as dis- 
cussed in Sect. 13.31 

The nuclear region interestingly falls into the LINER 
re gime. This contrast s with the optical classification as Syl.8 
by l Osterbrock ( 198ll) or from analysis of the SDSS spectrum 
(Fig. [TsT l. This discrepancy may be due to aperture effects, 
since the SDSS spectrum has been measured in a 3"diameter 
fiber covering the total FOV of our NIR observations. The total 
narrow-line flux of the H-recombination lines is dominated by 
the extra-nuclear regions. The LINER character might there- 
fore be hidden in the large-apertu re optical spectr um. Such ef- 
fects are also discussed by Kauffm ann et al.l (|2003|) . This can be 
tested in the present case by computing the flux ratio of the line 
fluxes integrated over the SINFONI FOV (~ 1.3kpc x 1.3kpc). 
Indeed, the position in the diagnostics diagram is shifted to- 
wards values typically found for Seyfert galaxies (Fig. fT4l i. 
This demonstrates the importance of high angular resolution 
in studying the nature of activity in galaxies. It is furthermore 
interesting that the off-nuclear emission regions populate a re- 
gion in the diagnostic diagram where AGN are usually found. 
This indicates that in these regions both shock and photoion- 
ization are important. 

Regions 2, 3, and 4 have line ratios typical of AGN, al- 
though the [FeII]/Pa/3 ratio is close to starburst values. Region 

5 has a 1 -OS (1)/Bry ratio among those of starbursts, although 
the [FeII]/Pa^ ratio is unconstrained because of the lack of 
J band spectral information. The overlay of the Pag with the 

6 cm radio map (Fig. |3] (d)) from lUlvestad & WilsonI (Il984 
shows radio emission in all five regions. The two-point spec- 
tral index a6cm/2icm ~ 1-3, assuming a power-law of the form 
Fy oc V is on the high end of values observed f or synchrotron 
radiation dUlvestad & Wilsonlll984l: iFilho et allEoOO) . Such a 
steep spectral index is comm only found in star-forming late- 
type galaxies ( ^Condonl ll992'). Together with the 1-0S(1)/Bry 
value, this indicates the starburst nature of the emission region 
5. As mentioned above, the three other off'-nuclear ratios popu- 
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Fig. 15. Rest frame SP SS spectrum (from SDSS DR4 



Adelman-McCarthyll2005h . The upper insets are close-ups of 
the HyS and Ha lines. The spectra are measured through a 
3"diameter fiber 



late the transition region from starburst to AGN in the diagram. 
At a projected distance of about 420 pc, regions 2 and 4 still 
show the strong influence of the (LINER) active nucleus. The 
jet-like extension of the 6 cm emission towards regions 2 and 
4 might also suggest the importance of outflows and associ- 
ated shocks that are able to excite [Fe II] and/or H2. Region 3 
at a projected distance of about 630 pc lies even closer to the 
starburtst regime. This is expected, since the nuclear influence 
should decrease at larger separations. 

Table 5. Emission line fluxes from SDSS spectrum (Fig.fTSl). 



Line 


flux 
[10^" Wm-2] 


FWHM 

[km s-'] 


[Oil] 3727 A 


1.9 = 


tO.l 


480 


Hy8 narrow 


0.80 = 


t0.02 


190 


[OIII] 5007A 


7.7 = 


tO.l 


400 


[O I] 6300A 


0.71 = 


t0.04 


365 


Ka narrow 


4.9 = 


tO.3 


220 


Ha broad 


8.1 = 


t04 


3170 


[Nil] 6585A 


4.5 = 


tO.3 


300 



Larkin et al 



1998 *) 



investi- 



In their study of LINERs, 
gated several multi-wavelength correlations. They found that 
LINERs and Seyferts do not follow the correlations between 
the IRAS 25 um to 60 lum inde x and Ha/Bry and [O I]/Ha found 
bv lMouri & Taniguchil(ll992h for starburst galaxies. The IRAS 



hot (100-150 K) dust. Seyfert galaxies typically show a shallow 
index. Therefore thi s index can be used as a Seyfert identifier 
jdeGriip etaljl992l) . LINERs show a wide spread in the IRAS 
index, whereas t hey only show little variation in the Ha/Bry and 
[Ol]/Ua ratios jLarkin et alJll998h . The values for Mrk 609 
(q'(25 : 60) = -1.9; [OI]/Ha * 0.15; Hz/Bry ^ 0.55 for the to- 
tal FOV or Hi/Bry * 1 .23 for the nuclear region) fall between 
the starburst correlation and the region occupied by LINERs. 
Note that the optical line-ratio is derived from the 3"diameter 
aperture SDSS spectrum, whereas the NIR line-ratio is calcu- 
lated for region 1 and for the total FOV. Nevertheless, a trend 
towards LINER values can be recognized at both wavelengths. 
Besides the aperture effects mentioned ab ove, considerable 
variability in the e mission has been detected ( Rudv et al]|l988 ; 
IPappa et aL 2002h. which can be com bined with the duty cycle 
hypothesis of lEracleous et alj ( 119951) . The authors propose re- 
curring accretion events that enhance the non-stellar continuum 
and generate Seyfert characteristics like broad emission lines. 
After the accretion event, the ionizing flux drops and the high- 
ionization states weaken, whereas the lower-ionization lines 
persist much longer, because of longer decay and light cross- 
ing times compared to the broad-line region. In this respect, 
Mrk 609 might approach a quiescent activity state in which 
the relative contribution of lower-ionization and shock-driven 
emission is enhanced with respect to photoionization. The star- 
burst, since farther away from the nucleus, remains constant. 
Unfortunately, no multi-epoch/multi-wavelength observations 
of Mrk 609 exist that could substantiate this scenario; and due 
to the calibration problems we cannot cornpare our l ine fluxes 
with the NIR measurements of Goodrich d 19901) and lRix et al.l 
(1990). 

The panchromatic properties of Mrk 609 reveal its com- 
posite nature. In Fig. [16] we plot the spectral energy distri- 
bution (SED) and compare it to averag e SEDs of starburst , 
LINER, and spiral galaxies (taken from ISchmitt et al1ll997h . 
Also shown is th e average SED of radio-quiet quasars from 



Elvis et al.l ( 119941) . It can be seen that the X- ray flux is nicely 



represented by an AGN-like contribution. Papp a et al.l (120021) 
estimate from the FIR luminosity that a starburst can only con- 
tribute up to 50% of the X-ray luminosity. But the AGN itself 
cannot explain the longer wavelength fluxes, which appear to 
be a combination of the three other contributions. Only star- 
bursts can produce the strong FIR fluxes. The LINER SED also 
fits nicely into this picture. Note that Mrk 609 is missing a sig- 
nificant big blue bump, which is also absent in low-lumino sity 
AGN and might be e xplained by low ac cretion-rate models (IHo 



19991). Furthermore. IWard et al.l ( 119881) found a strong correla 



tion between the broad Hor and the 2-10 keV X-ray luminosity 
as a measure for t he ionizing strength of the active nucleus. 



25 fim and 60 /im fluxes are sensitive to warm (30-50 K) and 



Pappa et al.l (120021) find, that the Ha emission from the large 
aperture spectrum is underluminous by a factor of about 40 and 
suggest a deficit of ionizing photons. This is agrees with the 
missing UV-bump in the SED. Another possibility is the vari- 
ability of the nonstellar continuum over the 15 years between 
optical and X-ray emission. 

The nuclear continuum and emission-line morphology is 
not symmetric and possibly resembles that of a weak nu- 
clear bar (Fig. |2d), |5a)). The emission peaks of the Paa 
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Fig. 16. Spectral energ y distribution of Mrk 609 cons tructed 
from literature values dRudy et al.lll988t jPollettaet anll996) . 
The yalues are compared w ith average SEDs f rom normal 
galaxies, sta rbursts, LINERs dSchmitt et alJll997b . and radio- 
quiet QSOs OEMs et al.lll994l) . 



map are aligned with the major and minor axes of the bar 
In addition, the slight elongation of the 6 cm radio emis- 
sion towards regions 2 and 4 (Fig. Od)) could be interpreted 
as a jet, which might in pa rt be responsible for the heat- 
ing of regions 2 and 4 (e.g. Ivan Breugel & Deyl Il993b . The 
interpretation is still not clear from our data. On the other 
hand, the clumpy structure of the ISM seen in the recom- 
bination line maps might also resemble a ring-like structure 
or parts of a nuclear spiral. Such str uctures have been un- 
covered in both the observation (e.g. lE rwin&SpaA3l2002t 
Martini et al. 2003; Emsellem et al. 2003 ) and sirnu lation (e.g. 
Englmaier & Shlosmanl i2000t iMacieiewskil 120041) of several 
quiescent and active galaxies. Torques associated with the gas 
response to the stellar-bar potential allow considerable amounts 
of gas to flow close r to nuclear regions (~ y r ' on time 
scales of ~ 10** vr: lEnglmaier & Shlosmanll2004l) where vis- 
cous processes take over the dis sipation of angular momen- 
tum (iDuschl & Strittmatter 12004 '). From the symmetrical ap- 
pearance with respect to the bar's minor axis, we might ex- 
pect a northeastern counterpart for region 3, which is missing, 
however, in the hydrogen emission-line maps. Nevertheless, the 
reddening map exhibits some additional reddening at the po- 
sition of the proposed counterpart. The velocity maps of the 
hydrogen recombination line do not show any dynamical be- 
havior, which is expected due to the face on view of the galaxy 
and the limited spectral resolution. Therefore, we are currently 
not in a position to favor any of the above models. 



4. Summary and conclusions 

Mrk 609 was chosen from a sample of AO-suitable optical 
counterparts of lu minous ROSAT X-ray AG N at low to inter- 
mediate redshifts dZuther et al.l |2004 l2005b . It is one of the 



lowest redshift objects in the sample and is therefore best- 
suited for initial integral field observations. 

We have presented first results on the circumnuclear struc- 
ture of the starburst/Seyfert composite Mrk 609 in the NIR us- 
ing imaging spectroscopy with SINFONI. The morphology is 
complex, and the continuum images reveal a bar-like structure. 
The distribution of hydrogen recombination emission (Paa, 
Pa/J, Bry) is clumpy and peaks at the tip where the potential 
bar meets the spiral arms and in regions along the minor axis. 
Bars can account for the angular momentum transfer necessary 
for fueling nuclear activity. The 6 cm VLA emission is also ex- 
tended in the latter direction. Whether this emission is associ- 
ated with a jet or with resonances in the bar potential, is unclear. 
The presence of nuclear broad Paor and [Si VI] are clear indi- 
cators of the accretion of matter onto a nuclear super-massive 
black hole. 

The analysis of molecular hydrogen and [Fe II] emission in- 
dicates the importance of shock heating, although X-ray heat- 
ing by the nucleus and non-thermal contributions is possi- 
ble. The distribution of molecular hydrogen follows the con- 
tinuum shape, while that of [Fell] is aligned with the minor 
axis of the continuum and with the radio contours, as well as 
with the H-recombination line peaks. The nucleus itself shows 
signs of LINER activity, which can be recognized by its high 
[Fe II]/Pa/J and H2/Br7 values. Our integral-field data clearly 
resolve the nuclear and starburst activity in the central kilo par- 
sec. Extinction appears to play no crucial role in this region, 
since the H-recombination line ratios are consistent with unred- 
dened case-B values. This is supported further by the strong 
Lya emission. However, already small amounts of extinction 
have considerable effects at visible wavelengths, as can be seen 
in the flocculent morphology of Fig. [T] Continuum reddening 
in the NIR, on the other hand, can already be caused by small 
amounts of warm/hot dust. The reddening of the continuum 
found toward the nucleus and the other regions is typical of 
AGN and star-forming regions where the dust reprocesses UV 
radiation from hot stars or the active nucleus. 

We found large amounts of cold molecular gas, which pro- 
vides the fuel for the star-forming activity. M-type giants are 
the dominating stellar population in the NIR, which indicates 
starburst ages of about 100 Myr (for a single burst). The nu- 
clear stellar absorption features are considerably diluted by the 
non-stellar continuum. 

Our LINER classification, together with the published vari- 
ability of the non-stellar NIR emission, might be explained by 
the duty-cycle hypothesis in which short-lived accretion events 
occur periodically and lead to the appearance of Seyfert fea- 
tures in the high state and of low-ionization (shock driven) 
emission features in the low state. 

The X-ray emission can account for some of the H2 ex- 
citation. This might also explain the higher H2/Bry ratios for 
the off-nuclear regions than expected for starbursts. The [Fe II] 
emission, however, does not seem to be excited by X-rays. 
Shock excitation along the southeast/northwest axis might be 
caused partly by a radio jet impinging into the ISM. The main 
driver for both species therefore appears to be shock excitation. 

Our results show that high spatial resolution is vital 
for dissecting the nuclear and starburst activity in AGN. 
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Morphological peculiarities, such as nuclear bars and nuclear 
starburst rings, are best detected at NIR wavelength, since the 
mass dominating stellar populations have their emission max- 
imum in the NIR, and dust extinction is smaller than at visi- 
ble wavelengths. In conjunction with high spectral resolution, 
it is possible to trace the stellar and gaseous nuclear kinemat- 
ics in order to find important constraints for dynamical mod- 
els describing the infall of matter into the nuclear regions (e.g. 
Fathi et al.luOOSb . This is a prerequisite for the study of com- 
posite systems in which nuclear star formation and Seyfert ac- 
tivity occur at similar levels of intensity. 
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o,K o,r 0.15 0.2 0.2S 0.3 0.3^ 4; S ^ M m "i^ 

Fig. 2. H + K emission line maps of recombination lines. They were created with QFITSVIEW by summing the flux over the 
emission Une and subtracting the median value of the neighboring continuum (a-c). All images were smoothed with a 3- pixel 
Gaussian. The values of the color bars are in units of 10"'' W m"^. (a) Narrow Paa; contours are calculated relative to the peak 
flux (90, 70, 50, 30%). (b) Bry overlayed with Paa contours, (c) He I 2.059 ;um overlayed with Paa contours, (d) Continuum 
image in a region without emission lines around 2.2 /zm. The panels are 3" x 3" each. North is up and east is left. Regions are 
indicated in panel (a) from which we extracted cummulative spectra. Regions 2 and 4 have a projected distance of about 420 pc, 
region 3 « 630 pc, and region 5 « 1.3 kpc to the center. 
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Fig. 3. H + K emission line maps (using the same procedure as for Fig.|2]i of: (a) l-OS(l) (contours: 90, 70, 50, 30, and 10% of 
peak flux), (b) [Fell] 1.664 fim (contours: 90, 70, 50, 30, 20, 10, and 5% of peak flux), and (c) [Si VI] (contours: 90, 70, 50, 30, 
and 10% of peak flux). The colorbar values are in units of 10 '^ W nT^. Panel (d) di splays an overlay of the Paa map from Fig. 
12 a) with the 6 cm VLA map in (magenta) contours (from luivestad & Wilsonl[l984t) . The blue ellipse in the lower right corner 
indicates the 6 cm beam size. The panels are 3" x 3" each. North is up and east is left. 
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Fig. 4. y-band continuum and emission line maps of: (a) 1 .25 jim continuum (contours: 90, 70, 50, 30, 20, and 10% of peak flux), 
(b) He I 1.083 fim (contours: 90, 70, 50, 30, 20, and 10% of peak flux), (c) Fell 1.257 fim (contours: 90, 70, 50, 30, and 15% of 
peak flux), and (d) Pa/3 overlayed with Paa in contours (see Fig. 12 a)). Color bar values are in units of lO"'** W m"^. The panels 
are 3x3 arcsec^ each. North is up and east is left. 



